System Size, Energy, Pseudo-rapidity, and Centrality Dependence 

of Elliptic Flow 
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Lhis paper presents measurements of the elliptic flow of charged particles as a function of pseu- 
dorapidity and centrality from Cu-Cu collisions at 62.4 and 200 GeV using the PHOBOS detector 
at the Relativistic Heavy Ion Collider (RHIC). The elliptic flow in Cu-Cu collisions is found to be 
significant even for the most central events. For comparison with the Au-Au results, it is found that 
the detailed way in which the collision geometry (eccentricity) is estimated is of critical importance 
when scaling out system-size effects. A new form of eccentricity, called the participant eccentricity, 
is introduced which yields a scaled elliptic flow in the Cu-Cu system that has the same relative 
magnitude and qualitative features as that in the Au-Au system. 



The characterization of the collective flow of produced 
particles by their azimuthal anisotropy has proven to 
be one of the more fruitful probes of the dynamics of 
heavy ion collisions at the Relativistic Heavy Ion Collider 
(RHIC). Flow is sensitive to the early stages of the colli- 
sion and so the study of flow affords unique insights into 
the properties of the hot, dense matter that is produced, 
including information about the degree of thermalization 
and its equation of state 0. In particular, flow measure- 
ments as a function of pseudorapidity have motivated the 
development of three-dimensional hydrodynamic models 
of relativistic heavy ion collisions [2j and provide infor- 
mation crucial in constraining these models and others 
that seek to better understand what role the longitudinal 
dimension plays in the collision 0, 0] . By examining the 
successes and failures of hydrodynamic models of heavy 
ion collisions, a better understanding of the dynamics of 
these collisions is gained. 

Elliptic flow has been studied extensively in Au-Au col- 
lisions at RHIC as a function of pseudorapidity ^.cen- 
trality, transverse momentum, and energy H, IE S 13 • 
The large pseudorapidity coverage of the PHOBOS de- 
tector makes it ideally suited for probing the longitudinal 
structure of the collision, the dynamics of which have 
only recently begun to be understood away from midra- 
pidity JlOj . This work presents new results on data taken 



by the PHOBOS experiment at RHIC showing a detailed 
comparison of differential measurements of elliptic flow 
in Cu-Cu and Au-Au collisions at ^/s^ = 62.4 and 200 
GeV. 

The strength of the elliptic flow, V2, is given by the co- 
efficient of the second harmonic in the Fourier expansion 
of the particle azimuthal angle distribution relative to the 
reaction plane, For this analysis ^2, an estimate of 
\l/ri, was used as in Ref. [llj . 

The PHOBOS detector is comprised of silicon pad de- 
tectors for tracking, vertex detection, and multiplicity 
measurements. Details of the setup and the layout of the 
silicon sensors can be found elsewhere |l| . Key elements 
of the detector used in this analysis include the first six 
layers of both silicon spectrometer arms, the silicon ver- 
tex detector (VTX) , the silicon octagonal multiplicity de- 
tector (OCT), three annular silicon multiplicity detectors 
to either side of the collision point, and two sets of scin- 
tillating paddle counters for centrality determination. 

Monte Carlo simulations of the detector performance 
were based on the HIJING event generator [l 7 | and the 
GEANT 3.211 [lj] simulation package, folding in the sig- 
nal response for scintillator counters and silicon sensors. 

The data shown here were taken with the PHOBOS 
detector at RHIC during the years 2001-2005. The Au- 
Au data are published in previous papers describing work 
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FIG. 1: V2 vs. 77 for Cu-Cu collisions at ^fs^ = 62.4 and 
200 GeV using the hit-based analysis. The boxes show the 
90% C.L. systematic errors and the bars represent the 1-er 
statistical errors. Previously published 200 GeV Au-Au data 
(without error bars) is shown for comparison. 



at ^/s^ = 19.6, 62.4, 130, and 200 GeV [HE El. The 
Cu-Cu data at ^/s NN = 62.4 and 200 GeV presented here 
were analyzed in a similar fashion, using the hit-based 
method that utilized hits in the VTX, OCT and ring 
sub-detectors to measure flow over a wide range in pseu- 
dorapidity (\r/\ < 5.4) and a track-based method that 
made use of tracks in the spectrometer arms and had a 
smaller pseudorapidity coverage (0.0 < 77 < 1.0). For 
details on the hit-based and track-based methods, see 
Refs. |5J and [6J, respectively. 

The event-by-event collision vertex was determined us- 
ing the intersection of tracks identified in the spectrome- 
ter and extrapolated back to a common point. The flow 
analysis was based on the anisotropy of the azimuthal 
distribution of charged particles traversing the detector. 
At the points where charged tracks passed through an 
active silicon detector, energy was deposited in the form 
of ionization. A pad where energy was deposited is said 
to be a "hit" . This analysis is based on the "sub-event" 
technique wherein one studies the correlation of hits in 
one part of the detector with the event plane angle as de- 
termined by hits in a different part of the detector [lT| . 
As described in the earlier work 0, 0, 0, 0] , corrections 
are applied to account for signal dilution due to detector 
occupancy and adjustments are made to create an ap- 
propriately symmetric acceptance for the analysis. The 
sub-event regions used in the event plane calculation were 
0.1 < |ry| < 3.0 for both 62.4 and 200 GeV. The event 
plane resolution was calculated separately for each cen- 
trality bin. The resolution correction factor ranged from 
2 to 3 on average, with the larger correction necessary 
at 62.4 GeV. For the determination of V2 in the positive 
(negative) 77 region of the detector, the sub-event on the 
opposite side of midrapidity was used to evaluate \&2- 



FIG. 2: V2 vs. ybcam for Cu-Cu collisions at ^/s^ = 62.4 
and 200 GeV from the hit-based analysis. Only 1-cr statistical 
errors are shown. 



Monte Carlo simulations showed a residual suppression 
of the flow signal dominated by background particles car- 
rying no flow information and the loss of sensitivity due 
to the hit map symmetrization and the occupancy correc- 
tion algorithm. As in our earlier work with the hit-based 
technique, this suppression was corrected using simulated 
data by comparing the output resolution corrected flow 
signal to the input flow signal for many samples of simu- 
lated data with different shapes and magnitudes of input 
flow signal. 

Numerous sources of systematic error were investi- 
gated, including effects due to the hit definition, hit merg- 
ing, sub-event definition, knowledge of the beam orbit 
relative to the detector, hole filling procedure, vertex- 
ing algorithm, and suppression correction determination. 
The effect of these sources depended both on 77 and cen- 
trality. In general, the systematic error arising from each 
source was determined by varying that specific aspect of 
the analysis (or several aspects in concert) within rea- 
sonable limits and quantifying the change in the final V2 
result as a function of 77 and centrality. The individ- 
ual contributions were added in quadrature to derive the 
90% confidence level error shown in the results presented 
here. The systematic uncertainty was dominated by the 
suppression correction determination. 

Figure ^ shows the elliptic flow signal as a function 
of pseudorapidity in Cu-Cu collisions at ^/s NN = 62.4 
and 200 GeV for the 40% most central events. The 
resemblance to published Au-Au results (also shown 
in Figure ^| is striking. The Cu-Cu V2 displays a sim- 
ilar shape in pseudorapidity to that of Au-Au, with a 
magnitude at midrapidity only 10-20% lower than that 
seen in Au-Au, increasing to ^40% at large \rj\. The 
strength of the Cu-Cu V2 signal is surprising in light of 
expectations that the smaller system size would result in 
a much smaller flow signal ^(| . 
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FIG. 3: V2 vs. Npart for Cu-Cu collisions at ^/s NN = 62.4 and 
200 GeV. The boxes show the 90% C.L. systematic errors and 
the lines represent the l-cr statistical errors. The results from 
two analysis methods are shown, as discussed in the text. V2 
is shown for \r]\ < 1 and < r\ < 1 for the hit-based and 
track-based methods, respectively. 



The Cu-Cu V2 also exhibits extended longitudinal scal- 
ing, as shown in Figure [2 and as already seen in Au-Au 
collisions for elliptic flow 0] and directed flow 01 an d 
for charged particle multiplicity 0, 0, |2(| . The agree- 
ment between the two energies in — j/beam implies that, 
as with Au-Au, the elliptic flow is largely independent of 
energy when viewed (effectively) in the rest frame of one 
of the colliding nuclei. 

The centrality dependence of the elliptic flow measured 
in Cu-Cu is presented in Figure^ where V2 is plotted as a 
function of the number of participating nucleons, N part . 
Both hit-based and track-based analyses were used for 
the 200 GeV data, and the results of the two methods 
agree quite well within errors. 

A substantial flow signal is measured in Cu-Cu at both 
energies for even the most central events. This is quite 
surprising, as the initial spatial anisotropy gives rise to 
a momentum space anisotropy which, in turn, produces 
the flow |2l|. ft is expected therefore that V2 should ap- 
proach zero as the collisions become more central, as it 
does for Au-Au |6( . The persistent and non-trivial elliptic 
flow signal seen in the most central events implies that 
something beyond the expected nuclear collision geome- 
try may be responsible for driving the flow signal. 

To explore this question further, it is useful to compare 
directly the elliptic flow signal across different colliding 
species, i.e., make a direct comparison between the flow 
seen in Cu-Cu and Au-Au collisions. To do this, it is 
important to scale out the difference in the initial geo- 
metric asymmetry of the collision, i.e., the eccentricity 
of the collision. This is crucial since in each selected 
centrality range the average eccentricity depends on the 
size of the colliding species. 



FIG. 4: The average eccentricity defined in two ways, ((e s td) 
and (epart}), as described in the text, vs. N part for simulated 
Cu-Cu and Au-Au collisions at ^/s NN = 200 GeV. The bands 
show the 90% C.L. systematic errors. 

Typically, the eccentricity is defined by relating the 
impact parameter of the collision in a Glauber model sim- 
ulation to the eccentricity calculated assuming the minor 
axis of the overlap ellipse to be along the impact param- 
eter vector. Thus, if the ir-axis is defined to be along 
the impact parameter vector and the y-axis perpendic- 
ular to that in the transverse plane, the eccentricity is 
determined bv plE^f 



where a x and a y are the RMS widths of the participant 
nucleon distributions projected on the x- and y-axes, re- 
spectively. Let us call the eccentricity determined in this 
fashion e s td- 

The relation of the eccentricity to the centrality de- 
pends on the details of the eccentricity definition used in 
the Glauber model simulation. The definition most com- 
monly used is presented above. Implicit in this choice is 
a physics bias about the relevant asymmetry that drives 
the flow signal. It is important to consider other possi- 
bilities. In particular, a natural choice to consider is the 
geometry of the participant nucleons themselves. 

In a large system, the nuclear geometry and the par- 
ticipant geometry largely coincide. For small systems or 
small transverse overlap regions, however, fluctuations in 
the nucleon positions in Glauber model calculations, as 
described below, frequently create a situation where the 
minor axis of the ellipse in the transverse plane formed 
by the participating nucleons is not along the impact pa- 
rameter vector. One way to address this issue is to make 
a principal axis transformation, rotating the x- and y- 
axes used in the eccentricity definition in the transverse 
plane in such a way that a x is minimized. Let us call the 
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FIG. 5: (a) v 2 (unsealed) vs. N part , (b) v 2 /(e s td) vs. N part , 
and (c) V2/(£ P art) vs. Npart, for Cu-Cu and Au-Au collisions 
at ^/s NN = 62.4 and 200 GeV. 1-a statistical error bars are 
shown. V2 is shown in \r/\ < 1 for the hit-based method. 



eccentricity determined in this fashion the participant 
eccentricity, e par t, and the plane specified by the beam 
axis and the ir-axis in the participant frame ^participant- 
In terms of the original x- and y-axes (in fact, any pair 
of perpendicular transverse axes), 



-part — 



(2) 



In this formula, a xy — (xy) — (x) (y) . The average values 
of e s td and e P art are quite similar for all but the most 
peripheral interactions for large species, as is shown in 
Figure Q] for Au-Au. For smaller systems such as Cu-Cu, 



however, fluctuations in the nucleon positions in Glauber 
model calculations (described below) become quite im- 
portant for all centralities and the average eccentricity 
can vary significantly depending on how it is calculated. 
This is also illustrated in Figure 0] 

The effects of finite number fluctuations on elliptic 
flow have been studied for large collision systems with 
Monte-Carlo simulations [24l l2a |26| and were found to 
be small. However, in Cu-Cu collisions these fluctua- 
tions are larger and could have a significant impact on 
the elliptic flow. The participant eccentricity allows these 
fluctuating configurations to be considered seriously on 
an event-by-event basis. 

The Glauber model used for the calculation of these 
eccentricities is a Monte Carlo toy model which builds nu- 
clei by randomly placing nucleons according to a Woods- 
Saxon distribution. Excluded volume effects were incor- 
porated into the model, requiring a minimum center- 
to-center nucleon separation of 0.4 fm, to agree with 
fflJING A number of sources of systematic er- 

ror were studied, including nuclear radius, nuclear skin 
depth, nucleon- nucleon inelastic cross-section ctnn, and 
minimum nucleon separation. The systematic error con- 
tributed by each source was determined by varying that 
specific parameter in the analysis within reasonable lim- 
its and quantifying the change in the final eccentricity 
result as a function of centrality. The individual contri- 
butions were added in quadrature to determine the 90% 
confidence level errors shown in Figure 0] 

The crucial importance of the definition of eccentric- 
ity in comparing Cu-Cu and Au-Au results can be seen 
in Figure \5\ where comparisons are made between Cu- 
Cu and Au-Au data using the eccentricity-scaled elliptic 
flow. In Fig. 03, V2/e s td increases rapidly in Cu-Cu as 
the events become more central, and is generally larger 
than that of Au-Au. One might then conclude from this 
that either the smaller Cu-Cu system produces V2 much 
more efficiently than the larger Au-Au system or that 
e s td may not be the appropriate quantity for describ- 
ing the initial geometry of the collision. Consider then 
Fig. |5j;, in which v 2 /£p art is shown to be very similar 
for both Cu-Cu and Au-Au, even appearing to lie on the 
same curve. Given the qualitative and quantitative sim- 
ilarities between the results in the two systems already 
shown, it is not unreasonable to expect both systems 
to have a similar eccentricity-scaled elliptic flow, as in 
Fig-Efc- Therefore, it seems likely that e part as discussed 
— or a rather similar quantity, such 

is the relevant eccentricity for the 



here and in Ref. 

as J(e* t ) & 



~partl 

azimuthal anisotropy. 

In summary, the results presented here show a measur- 
able and significant elliptic flow signal in Cu-Cu collisions 
at 62.4 and 200 GeV. These data show that qualitative 
features attributed to collective effects in Au-Au persist 
down to the relatively small numbers of participants seen 
in the Cu-Cu collision and are of comparable magnitude. 
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The essential role of the choice of collision geometry in 
comparing flow across nuclear species is clearly demon- 
strated. If it is assumed that the flow is independent 
of species for a given collision geometry, the data shown 
here strongly suggest that it is the participant eccentric- 
ity, not the nuclear eccentricity, that is responsible for 
elliptic flow. These results also imply that ^2 is an es- 
timate of ^participant, which may be the relevant angle, 
particularly for systems with smaller numbers of partic- 
ipants. This, in turn, may provide information on the 



nature of the matter driving the flow. 
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